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ABSTRACT 

We present the first results obtained from tiie identification of ~ 30 000 sources in the Spitzer/24: fim 
observations of the COSMOS field at S24f_im ^ 80 /iJy. Using accurate photometric redshifts (ctz ~ 0.12 
at z ^ 2 for 24 /im sources with i+ < 25 mag AB) and simple extrapolations of the number counts 
at faint fluxes we resolve with unprecedented detail the build-up of the mid-infrared background 
across cosmic ages. We find that ~ 50% and ~ 80% of the 24 /im background intensity originate from 
galaxies at 2; < 1 and z <2 respectively, supporting the scenario where highly obscured sources at 
very high redshifts (z > 2) contribute only marginally to the Cosmic Infrared Background. Assuming 
flux-limited selections at optical wavelengths, we also find that the fraction of i+-band sources with 
24 fim detection strongly increases up to z ~ 2 as a consequence of the rapid evolution that star- 
forming galaxies have undergone with lookback time. Nonetheless this rising trend shows a clear 
break at z~1.3, probably due to /c-correction effects implied by the complexity of spectral energy 
distributions in the mid-infrarcd. Finally, we compare our results with the predictions from different 
models of galaxy formation. We note that semi-analytical formalisms currently fail to reproduce 
the redshift distributions observed at 24 fim. Furthermore the simulated galaxies at «S'24^m > 80 /iJy 
exhibit R — K colors much bluer than observed and the predicted i^-band fluxes are systematically 
underestimated at z>0.5. Unless these discrepancies mainly result from an incorrect treatment of 
extinction in the models they may reflect an underestimate of the predicted density of high redshift 
massive sources with strong on-going star formation, which would point to more fundamental processes 
and/or parameters (e.g., Initial Mass Function, critical density to form stars, feedback, ...) that are 
still not fully controlled in the simulations. The most recent backward evolution scenarios reproduce 
reasonably well the flux/redshift distribution of 24 /im sources up to 2;~3, although none of them is 
able to exactly match our results at all redshifts. 

Subject headings: galaxies; high-redshift — infrared: galaxies — cosmology: observations 



1. INTRODUCTION 

Over the last decade, extragalactic infrared (IR) sur- 
veys have become a key component in our quest to un- 
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derstand galaxy formation. The reprocessing of stel- 
lar and nuclear radiations by dust and the result- 
ing emission of this energy at thermal IR wavelengths 
(8 /im < A < 1000 /im) plays a major role in shaping the 
panchromatic appearance of star-forming galaxies and 
active galactic nuclei (e.g., Silva et al. 1998; Gordon et al. 
2000; Nenkova et al. 2002; Dopita et al. 2005; Marshall 
et al. 2007; Siebenmorgen & Kriigel 2007; da Cunha et al. 
2008). The direct characterization of the properties of 
these sources at IR wavelengths is therefore crucial for 
reliably estimating the fraction of their bolometric en- 
ergy output absorbed by dust and converted into lower 
energy photons. It also provides unique constraints on 
the nature of the different physical processes driving their 
on-going activity (e.g, de Grijp et al. 1985; Mazzarella 
et al. 1991; Genzel et al. 1998; Dale et al. 2006; Smith 
et al. 2007; Spoon et al. 2007). 

The discovery of a large number of distant sources ra- 
diating copious amounts of energy in the IR (e.g.. Small 
et al. 1997; Barger et al. 1998; Aussel et al. 1999; El- 
baz et al. 1999; Blain et al. 2002; Papovich et al. 2004; 
Dole et al. 2004; Daddi et al. 2005) revealed that critical 
phases of the cosmic star formation and nuclear accretion 
history happened within dust-embedded environments 
hidden behind large amounts of absorbing material. A 
significant fraction of the growth of structures over cos- 
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mic ages thus occurred in galaxies affected by substantial 
reddening and which characterization at shorter wave- 
lengths (e.g., UV, optical) can be severely biased if their 
extinction is not properly estimated. At large cosmo- 
logical distances, highly obscured sources can even be 
barely detectable below ~ 1 /zm and the observations in 
the IR or the millimeter are sometimes the only way to 
unveil the presence of on- going star formation and/or nu- 
clear activity in these galaxies (e.g., Hughes et al. 1998; 
Downes et al. 1999; Houck et al. 2005). 

The infrared spectral energy distributions (SED) of 
star-forming galaxies and active galactic nuclei (AGN) 
usually peak at ^ 80-200 /im. The most reliable esti- 
mate of the energy absorbed by dust in these sources 
would thus require observations in the far-IR or the sub- 
millimeter. Unfortunately though, data taken at such 
wavelengths have been so far hampered by source con- 
fusion as well as the very modest sensitivity of current 
detectors and the difficulty to access this frequency range 
from the ground. For instance, observations of the dis- 
tant Universe performed at 70 /xm and 160 ^m with the 
Spitzer Space Telescope (Werner et al. 2004) or in the 
submillimeter with ground-based cameras of bolometers 
have been so far limited to the identification of galaxies 
at the very bright end of the luminosity function (i.e., 
Lboi^ 10^^ L0), with biases due to dust temperature se- 
lection effects (e.g., Blain et al. 2004; Frayer et al. 2006; 
Pope et al. 2006; Borys et al. 2006; Le Floc'h et al. 2007). 
Detecting the more typical high redshift sources at the 
peak of their IR SED has not been achieved yet, and our 
knowledge of the far-IR luminosity function of galaxies 
in the distant Universe is still affected by substantial un- 
certainties. 

This strong limitation can be partially remedied with 
observations performed at mid-IR wavelengths (i.e., 
5 /im < A < 60 /im) . Although cameras operating in the 
mid-IR are not as efficient as optical detectors they have 
been far more sensitive than in the far-IR, allowing the 
probe of the IR luminosity function down to much fainter 
levels. The observed correlations between the mid-IR 
and the far-IR luminosities of star-forming galaxies al- 
low reasonable estimates of their total IR luminosity 
(LiR = Lg-ioooAim) with extrapolations from fluxes mea- 
sured in the mid-IR (Spinoglio et al. 1995; Chary & El- 
baz 2001; Bavouzet et al. 2008; Symeonidis et al. 2008). 
At such wavelengths, starbursts and AGNs exhibit also 
quite disctinct spectral signatures, providing invaluable 
clues to quantify the relative contribution of star-forming 
and nuclear activity occurring within individual sources. 

Observations of the deep Universe with the 24 /xm 
bandpass of the Multi-band Imaging and Photometer for 
Spitzer (MIPS, Ricke et al. 2004) and with the Infrared 
Spectrograph (IRS, Houck et al. 2004) have led to tremen- 
dous progress in our understanding of the bolometric 
properties of distant sources (e.g., Huang et al. 2005; 
Yan et al. 2005; Daddi et al. 2005; Papovich et al. 2006; 
Webb et al. 2006; Reddy et al. 2006; Menendez-Delmestre 
et al. 2007; Sajina et al. 2007; Lutz et al. 2008; Pope et al. 
2008). The Spitzer data provide tight constraints on the 
different mechanisms (merging, feedback, environmental 
effects) that regulate the activity of star formation and 
nuclear accretion in galaxies (e.g., Noeske et al. 2007; El- 
baz et al. 2007; Bridge et al. 2007; Buat et al. 2008; Bai 



et al. 2007; Marcillac et al. 2008). They also enable the 
identification of a large number of obscured AGNs that 
must have played a key role in driving the coeval growth 
of bulges and super-massive black holes but whose contri- 
bution is usually underestimated by surveys performed 
at other wavelengths (e.g., Houck et al. 2005; Polletta 
ct al. 2006; Alonso-Herrero et al. 2006; Dey et al. 2008; 
Fiore et al. 2008). 

In spite of this wealth of results recently achieved with 
Spitzer at mid-IR wavelengths, a coherent and compre- 
hensive picture of the deep IR Universe is still missing 
though. Given the variety of starburst and AGN SEDs 
in the mid-IR (e.g., Laurent et al. 2000; Brandl et al. 
2006; Wcedman et al. 2005; Wu et al. 2006; Armus ct al. 
2007), observations at these wavelengths are affected by 
quite complex fc-correction effects and their sensitivity 
decreases rapidly beyond z^l. Hence, wide and shal- 
low surveys are only sensitive to the very bright end of 
the mid-IR luminosity function, while the deeper obser- 
vations usually performed at the expense of the covered 
area suffer from limited statistics. Also, the identifica- 
tion of mid-IR selected sources at high redshift can be 
extremely difficult given their faintness at optical wave- 
lengths. As a result the exact role that dusty galaxies 
played with respect to the less luminous but much more 
numerous sources at the faint end of the luminosity func- 
tion is still under debate (but see Reddy et al. 2008). 

In the attempt to unify the different pictures of galaxy 
evolution recently achieved by optical and IR surveys we 
undertook deep 24 /xm imaging over the whole area of 
the COSMOS field using the MIPS instrument on-board 
Spitzer. The Cosmic Evolution Survey (COSMOS) is de- 
signed to probe the evolution of galaxies and AGNs up to 
z ^6 over a sky region large enough to address the role 
of environment and large scale structures (Scoville et al. 
2007b). It is based on deep multi- wavelength observa- 
tions performed over ^ 2 deg^ from the X-ray to radio 
wavelengths (Hasinger et al. 2007; Taniguchi et al. 2007; 
Sanders et al. 2007; Bertoldi et al. 2007; Schinnerer et al. 
2007), with additional support from high spatial resolu- 
tion imaging obtained with the Hubble Space Telescope 
(Scoville et al. 2007a; Kockemoer ct al. 2007). Inten- 
sive spectroscopic follow-up was also carried out with 
the Very Large Telescope (Lilly et al. 2007) while photo- 
metric redshifts with unprecedented accuracy were deter- 
mined for more than ~ 900 000 z+-band selected galax- 
ies in the field (Ilbert et al. 2009). The quality of the 
data sets available in COSMOS (and in particular the 
depth with respect to the covered area) is therefore very 
well suited for probing the multi-wavelength properties 
of galaxies with minimized cosmic variance and over a 
wide range of redshifts and luminosities. 

In this paper we present our very first analysis of 
the mid-IR selected galaxy population detected down 
to <S'24/im ~ 80 yuJy in COSMOS. Data are described in 
Sect. 2 and our 24 /xm source number counts are shown 
in Sect. 3. Using the photometric redshifts of Ilbert et al. 
(2009) and Salvato et al. (2009) we then determine the 
24 /im redshift distributions in Sect. 4 and we analyze 
the number counts as a function of cosmic time in Sect. 5 
to constrain the history and the build-up of the Cos- 
mic Infrared Background (CIB). In Sect. 6 we present 
comparisons between the observed 24 /xm source redshift 
distributions and the predictions from both phenomeno- 
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logical and semi-analytical scenarios of galaxy formation. 
The results arc discussed in Sect. 7, which allows us to 
quantify how well the contribution of dusty luminous 
sources as a general population is currently understood 
in models of galaxy evolution. Future papers will be ded- 
icated to the analysis of the multi- wavelength properties 
of the COSMOS mid-IR selected galaxies, their evolu- 
tion, their clustering and their connection with sources 
selected with other criteria. Throughout this paper we 
assume a ACDM cosmology with Hq = 70 km s~^ Mpc^^ , 
fi„ = 0.3 and Qx = 0.7'(Spergcl et al. 2003). Mag- 
nitudes are quoted in AB^** unless the use of the Vega 
system is explicitly stated (e.g.. Sect. 6.1). 

2. THE DATA 
2.1. The 24 ^irn observations of COSMOS 
2.1.1. Data acquisition and reduction 

The COSMOS field was observed at 24 /im with the 
MIPS instrument on-board Spitzer as part of two Gen- 
eral Observer programs (PI D. Sanders). The first obser- 
vations (G02, PID 20070) were carried out in January 
2006 with the MIPS medium scan mode and two scan 
passes per sky pixel. As described in more detail by 
Sanders et al. (2007) this shallow imaging covered a to- 
tal area of ~ 4 deg^ centered on COSMOS, with a median 
integration time of 80 s per pixel. This first program was 
also used to test in a small region of the field the feasi- 
bility to improve the 24 /xm sensitivity with longer inte- 
grations despite the COSMOS mid-IR background being 
~2x higher than in other typical "cosmological" fields 
(e.g., Hubble/ Chandra Deep Fields, Lockman Hole). Af- 
ter validating this additional step, deeper imaging was 
performed in 2007 over the whole COSMOS field as part 
of a second program using the MIPS slow scan mode 
(G03, PID 30143). These data yielded a second and in- 
dependent 24 /xm coverage over a total area of >3deg^. 
Across the nominal 2 deg^ of COSMOS the combination 
of the G02 and G03 observations results in a median 
integration time of ~ 3 360 s per sky pixel. Note that the 
MIPS scan mode provides simultaneous imaging at 24, 
70 and 160 /zm; the observations of the COSMOS field at 
70/160 /im are presented by Frayer et al. (submitted to 
AJ). 

All the data set was reduced following the procedure 
outlined in Sanders et al. (2007). For each observing 
campaign a background pattern derived for every scan 
mirror position of the instrument was scaled and sub- 
tracted from all individual frames. These background- 
siibtracted images were then interpolated to a common 
grid of pixels using the MOPEX package (Makovoz & 
Khan 2005). They were finally coadded with our own 
IDL routines using the median of all the pixels associ- 
ated with each sky position. We stress that the use of 
a median combination was necessary in order to reject 
moving sources present in the field at the time of the 
observations. A forthcoming paper (Aussel et al. in 
prep.) will describe in more detail the data reduction. 
It will also provide a list of the 24 /um-detected asteroids 
in COSMOS. 

2.1.2. 24l^m source extraction, photometry 
1" mag (AB) = -2.5 x logio [F,. (mjy)] +16.4. 



The 24 /im sources in our mosaic were detected us- 
ing the automatic procedure of the Sextractor software 
(Bertin & Arnouts 1996) and we measured their flux den- 
sities with the Point Spread Function (PSF) fitting tech- 
nique of the DAOPHOT package (Stetson 1987). Given 
the Full- Width at Half Maximum (FWHM) of the MIPS- 
24 /im PSF 6"), most of the sources are not resolved in 
our data. In deep and crowded images such as the ones 
typically obtained in the MIPS deep surveys, the PSF 
fitting thus provides more reliable fiux measurements 
than aperture photometry. Furthermore, DAOPHOT 
performs simultaneous fits to multiple objects, which en- 
sures a proper separation of blended sources especially 
when the fluxes are close to the confusion limit. No prior 
information from any other wavelength was used for the 
source detection itself. However, several PSF fitting iter- 
ations were necessary to achieve a clean extraction of the 
blended sources that were not properly identified by Sex- 
tractor in the input 24 /im image. In brief, these blended 
cases were located (with Sextractor again) as positive 
detections in the residual mosaic obtained after the first 
PSF fitting. Their positions were added to the initial 
list of 24 /im sources determined with Sextractor and an- 
other extraction of the initial COSMOS 24 /xm image was 
then performed. Three iterations of this process were 
carried out, until the residual mosaic was cleaned from 
any positive signal above the detection threshold chosen 
for Sextractor. The result of this PSF fitting is illus- 
trated in Figure 1. It shows a sub-region of our COSMOS 
24 yum mosaic along with the corresponding residual im- 
age obtained after subtracting all sources extracted with 
DAOPHOT. 

To run this PSF fitting we built an empirical PSF from 
the brightest point sources identified in the mosaic and 
we subsequently fit this PSF to every single object de- 
tected by Sextractor. DAOPHOT performs this fit by 
scaling the empirical PSF within a fixed circular aper- 
ture centered on the source and it provides the flux of 
the fitted PSF enclosed within this aperture. Therefore a 
correction must be applied to account for the spatially- 
extended wings of the PSF lying beyond the radius of the 
chosen aperture. This correction was determined using 
the MIPS-24 /im STiny Tim PSF models provided by the 
Spitzer Science Center. These models allow to accoimt 
for the fraction of energy lying at large radii where the 
signal to noise in our data is not large enough to constrain 
the exact profile of the Point Spread Function. Relying 
on such corrections we believe that our flux absolute cali- 
bration should be accurate within a few percent, which is 
consistent with the results obtained by Engelbracht et al. 
(2007) and Rieke et al. (2008) as part of the calibration 
of the MIPS 24 /im detector. 

Following the convention adopted by the Spitzer Sci- 
ence Center we assumed a stellar 10 000 K black-body 
spectrum as the reference underlying SED for the 24 /im 
flux density measurements. We tested this calibration 
by cross-correlating our source list with the stars of the 
Two Micron All Sky Survey (2MASS) catalog (Jarrett 
et al. 2000). Our comparison did not reveal any sys- 
tematic deviation from the relation Kg — [24] = (Vega, 
with [24] denoting the equivalent magnitude measured in 
the 24;um MIPS bandpass), which is expected for most 
of the stellar spectral types when no dust emission ex- 
cess is detected above the photospheric stellar emission 
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Fig. 1. — A 15'xl5' sub-region of the COSMOS field observed at 24 /im (left panel), along with the residual image obtained on the 
same field of view after subtracting sources with a PSF fitting technique (central panel). The "reconstructed" image derived with this PSF 
fitting is shown on the right panel. 



(e.g., Gorlova et al. 2006). This cross-correlation also al- 
lowed us to register the absolute astrometric calibration 
of our 24 /xm mosaic to the 2MASS astrometric system. 
This was performed by applying a systematic offset of 
6 Dec = — 0.3" to the astrometry derived from the Spitzer 
data products. 

The noise in the final 24 /xm mosaic was characterized 
using two different methods. First, we computed the dis- 
persion of a set of random flux measurements performed 
over blank field regions of the COSMOS 24 /im image 
within the same aperture as the one used for the PSF fit- 
ting. Similar to the photometry measurements described 
above a flux correction was applied to account for the ex- 
tent of the PSF lying outside of the aperture, which led to 
an equivalent sensitivity limit of la = 18 /xJy. Second, we 
also performed a number of simulations adding artificial 
point-like sources in blank field regions of the mosaic. 
We measured the fluxes of these sources with the PSF 
fitting of DAOPHOT and the distribution of the differ- 
ences between their measured and input fluxes was fitted 
with a gaussian function. Our best fit was obtained for a 
la standard deviation varying between 14 /xJy and 16 fxJy 
depending on the range of input fluxes. Using these sim- 
ulations we finally examined the robustness of the error 
measurements provided by our PSF fitting by comparing 
the DAOPHOT uncertainties obtained for the simulated 
point sources with the errors directly estimated from the 
differences between their measured and input fluxes. We 
applied a systematic scaling of 3.8 to all the DAOPHOT 
errors so as to get 68% of the distribution of these differ- 
ences within the DAOPHOT la uncertainties. 

The completeness of our 24 /xm source extraction was 
determined using the same point source simulations as 
the ones described in the previous paragraph. We in- 
serted within the mosaic a set of empirical PSFs that 
were scaled to provide a flux distribution representa- 
tive of the one characterizing the 24 /xm extragalactic 
source population. These PSFs were subsequently ex- 
tracted following the same method as the one employed 
to detect and to measure the flux of the sources already 
present in the 24 ^m image. We found that our source 
extraction is more than 90% complete above a 24 /xm flux 



'S'24^m ~ 80 /xJy, which should thus be considered as a safe 
and conservative flux limit for any analysis requiring an 
unbiased 24 /xm selection of sources in the COSMOS fleld. 
According to the simulations though, our PSF fltting is 
still reliable down to fluxes as faint as ^60 /iJy and in 
spite of a lower completeness (~ 75%) our source extrac- 
tion should suffer from a relatively small contamination 
by fake objects at these flux levels. 

Based on this characterization we defined our final 
COSMOS MlPS-24^m catalog as the population of 
24 /xm sources brighter than 60 /xJy lying within the area 
covered by the COSMOS photometric redshift catalog 
of Ilbert et al. (2009, see section 2.3) but outside of the 
masked regions defined from the Subaru optical observa- 
tions of the COSMOS field. These masked areas cover a 
surface as large as 0.26 dcg^ but the optical photometry 
in these regions is affected by larger uncertainties due to 
the presence of very bright or even saturated objects. In 
particular, the photometric redshifts in these areas are 
less accurate. Keeping these regions would thus compro- 
mise our analysis of the MlPS-selected galaxies given the 
large fraction of 24 /xm sources associated with faint op- 
tical counterparts. Our 24 /xm catalog represents a total 
of 39 413 sources covering an effective area of 1.68 deg^. 
In the remaining of this Section 2 we will describe the 
identification of the 24 /xm source counterparts down to 
'S'24/im = 60 /tJy. This characterization will supplement 
other studies carried out by our team using selections 
performed at other wavelengths. However, the analysis 
reported in some other sections of this current paper as- 
sume an unbiased 24 /tm selection and only sources above 
80 /tJy will be considered in this case. 

2.2. Identification of 24 /xm source counterparts 

To ensure the most reliable and most complete identi- 
fication of the 24 /im source counterparts, we first corre- 
lated the 24 /im data with the Ks band COSMOS cata- 
log of McCracken et al. (submitted to ApJ). As we will 
see later a non-negligible fraction of the 24 /xm sources at 
z > 1 are associated with extremely faint optical coun- 
terparts. A direct cross-correlation between the MIPS- 
24 /xm and the optical catalogs of COSMOS would thus 
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result in a lower rate of identifications. Furthermore, the 
larger uncertainties affecting the determination of source 
centroids at faint 24 /xm fluxes, the width of the 24 /xm 
PSF and the very high density of optical sources in our 
data (e.g., ~45 galaxies arcmin"^ for i^3<25mag in 
the Subaru COSMOS observations, Taniguchi et al. 2007; 
Capak et al. 2007) could lead to lots of wrong associa- 
tions if these identifications had to rely on a direct match 
between the 24 /im and optical source catalogs. On the 
contrary, the COSMOS Kg-h&nd catalog is deep enough 
{Kg ~ 23.7 mag at 5(t, McCracken et al.) to detect coun- 
terparts for the majority of the 24 /im selected sources 
down to our MIPS sensitivity limit. The density of 
Kg-hand sources is also substantially smaller 7 galax- 
ies arcmin"^ for Kg < 22 mag) than in the optical bands 
and the near-IR wavelengths are "closer" to the 24 /xm 
bandpass, which minimizes the risk of wrong associations 
between mid-IR and i^^-selected sources. Finally we 
note that these if^-band observations of COSMOS were 
executed under excellent seeing conditions (~0.7" at 
2.2 /im). Compared to the other COSMOS near-IR data 
taken with e.g., the IRAC instrument (FWHM ^ 1.6" at 
3.6/im, Sanders et al. 2007), the PSF in the X^-band 
image is much narrower and allows more robust identifi- 
cations in the case of sources blended at 24 /xm. 

We correlated the 24 /xm and ii'^-band catalogs with 
a matching radius of 2". We found that this distance 
was large enough to enable the identification of most of 
the 24 /im sources, taking into account the width of the 
24 /im PSF (FWHM ~ 6") and the uncertainty on the 
astrometric position of each single 24 /im source. This 
radius was nonetheless sufficiently small to minimize the 
fraction of multiple matches as well as the association be- 
tween objects randomly-aligned on the sky. Among the 
39 413 MIPS-24 /im sources of the catalog, we identified 
a single Kg-hand counterpart for 33 146 objects (84% of 
the sample) while two possible matches were obtained for 
4470 other 24 /im detections (~ 11.5%). For theses cases 
we kept the closest possible counterpart when its centroid 
was at least twice closer to the 24 /im source than the sec- 
ond possible match (2 320 objects) while additional pri- 
ors from the COSMOS IRAC-3.6/im observations were 
considered to assign the most likely counterpart to the 
other 24 /im detections. No association was foimd for 
the remaining 1797 sources (4.5% of the catalog). A 
cross-correlation between these unidentified sources and 
our 3.6 /im catalog using the same matching radius re- 
vealed IRAC counterparts for 437 objects, a rare pop- 
ulation that must be characterized by extremely steep 
SEDs between 2 /im and 24 /im. For the other sources a 
visual inspection of our data showed that some of them 
may have counterparts either fainter than the detection 
thresholds of our Kg /IRAC COSMOS catalogs or with a 
centroid located slightly beyond our 2" matching radius. 
Many of the unidentified objects have also a signal/noise 
below 5(7 at 24 /im. At the faintest 24 /im flux leves some 
could originate from a small contamination of our catalog 
by fake sources. It should however not significantly bias 
our subsequent results given the small fraction that they 
represent among the whole sample of 24 /im sources. In 
summary we identified secure Ks-h&nd counterparts for 
~90% of the COSMOS 24 /im catalog down to 60/iJy, 
while less robust identifications were still found for an- 
other 5.5% of the MIPS-24 /im source population. 



Finally we note that this cross-identification with near- 
IR data also revealed 53 pairs of 24 /im sources associated 
with the same ifg-band counterpart. A careful inspec- 
tion of the 24 /im mosaic revealed that they correspond 
to single objects incorrectly deblended by the PSF fit- 
ting performed in the MIPS data. These "double" com- 
ponents were merged together to form a unique 24 /im 
source, reducing the MlPS-selected catalog to 39 360 ob- 
jects. 



2.3. Photometric redshifts 

2.3.1. Input catalogs 

The distance of the 24 /im sources was characterized 
using the photometric redshifts derived by Ilbert et al. 

(2009) and Salvato et al. (2009) for the optically and 
X-ray selected sources of the COSMOS field. These pho- 
tometric redshifts were determined with deep photome- 
try (e.g., Capak et al. 2009, in preparation) performed 
in a total of 30 large, medium and narrow band filters 
over a wavelength range covering the far-UV at 1550A 
up to the mid-IR at 8.0 /im. This large amoimt of data 
and the high quality of the obsca'vations used in their 
work resulted in a substantial improvement of photomet- 
ric redshift accuracy with respect to the determination 
of redshifts performed in other fields. 

The catalog of Ilbert et al. (2009, hereafter 109) 
provides redshifts for 933 789 sources selected at 
i^B<26.5mag from the Subaru/Suprime-CAM obser- 
vations of COSMOS, with a dispersion as small as 
(Jaz/{i+z) = 0.012 for i^g<24mag and 2; < 1.25. This 
unprecedented accuracy in photometric redshift determi- 
nations was mostly achieved thanks to the high sampling 
of the observed spectral energy distribution of sources 
around the position of their Balmer break as well as a 
new method accounting for the contribution of emission 
lines to their measured flux densities. As shown by Il- 
bert et al. (2009) the imcertaintics affecting their redshift 
estimates depend primarily on the redshift and the ap- 
parent i^-band magnitude of sources (see their flgure 9), 
with errors naturally increasing for fainter and more dis- 
tant galaxies. However a comparison with faint spectro- 
scopic samples obtained in the COSMOS field (Lilly et al. 
2007) revealed a dispersion <7£^z/(i+z) as low as 0.06 for 
sources with 23 mag < < 25 mag at 1.5 < 2; < 3. These 
small uncertainties will be critical for our analysis of the 
MlPS-selected galaxy population given the large fraction 
of 24 /im sources associated with faint (i.e., i > 24 mag) 
optical counterparts at z>l. 

The redshifts provided by 109 were determined us- 
ing a library of star-forming and passive galaxy tem- 
plates and they are not adapted to galaxies dominated in 
their optical/near-IR continuum by the contribution of 
an active galactic nucleus (such as the type 1 AGNs and 
optically-luminous quasars). This issue must be treated 
with care given the non-negligible contribution of AGNs 
to IR-selected galaxy samples (e.g., Gcnzcl et al. 1998; 
Houck et al. 2005). The redshifts of the 24 /im sources 
also identified with an X-ray counterpart were therefore 
supplemented using the catalog of Salvato et al. (2009), 
which provides photometric redshifts for 1542 X-ray se- 
lected sources detected in the XMM-Newton observa- 
tions of COSMOS (Hasinger et al. 2007; Brusa et al. 
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2007). Salvato et al. (2009) used an optimized set of 
templates with hybrid combinations of AGN and non- 
active galaxy SEDs, which allowed them to derive pho- 
tometric redshifts with an accuracy comparable to that 
achieved by Ilbert et al. (2009). They also made an ex- 
tensive use of multi-epoch observations of COSMOS to 
correct for the time-variability of type 1 AGNs. The 
photometric redshifts of the full XMM-selected sample 
mostly lie at 0.5 < ^ < 1.5 but extend up to z ~ 3, with an 
accuracy as good as crAz/(i+z) = 0.015 for faint sources 
(ijB> 22.5 mag) at z~l (Salvato et al. 2009). 

We stress that the redshifts of 109 were also used in 
the case of the 24 /xm sources associated with obscured 
nuclei (such as sources characterized by a rising power- 
law continuum in the IRAC bands) but with no X-ray 
detection. Their emission at optical wavelengths is usu- 
ally dominated by the contribution of their host galaxy 
and the photometric redshifts of 109 should not be sig- 
nificantly less accurate than those determined for star- 
forming sources. This argument is actually reinforced by 
the rather good agreement that we observed betweeen 
the redshifts of Salvato et al. (2009) and those of Ilbert 
et al. (2009) for the type 2 AGNs common to both sam- 
ples. 

2.3.2. Photometric redshifts of 24 t^rn, sources 

To derive the photometric redshifts of the 24 /im 
sources we first correlated the list of their Ks band 
counterparts with the catalog of Ilbert et al. (2009) us- 
ing a matching radius of 1". We found a single opti- 
cal counterpart at i^g< 26.5 mag for 35 347 i^^-band 
sources (90% of the full 24 ^m source catalog). Double 
matches were obtained for 361 other Xg-band objects. 
In those cases, we noticed that the difference between 
the matching distances of the first and the second clos- 
est matches was large enough to consider the first one as 
the most likely optical counterpart. The catalog of 109 
is limited to ijg = 26.5 mag and no optical match could 
be found for the remaining 1855 JsTg-band sources (i.e., 
iX^ > 26.5 mag). 

From this match between the Ks and the optical bands 
we then identified 1 129 optical counterparts also associ- 
ated with a source detected by XMM. Among those, 
1111 are included in the analysis presented by Salvato 
et al. (2009) and we thus assigned to these objects the 
photometric redshifts provided in their catalog of X-ray 
sources. The identification of the optical counterparts 
associated with the remaining 18 XMM objects is con- 
sidered less secure (Brusa et al. 2007). For these objects 
we kept the original photometric redshift derived by 109. 

Finally we considered the list of the 1 797 MIPS-24 /im 
sources for which we did not find any iiTg-band coun- 
terpart (see Sect. 2.2). We correlated this list with the 
catalog of Ilbert et al. (2009) using a matching radius 
of 2", assuming that some of these sources could have 
counterparts detected in the very deep Subaru i^-band 
observations of COSMOS (Taniguchi et al. 2007) but not 
reported in our Ks band catalog. Among these 24 /im 
sources, 667 objects were identified with one single coun- 
terpart in the optical catalog. We note that some of 
them seem qto be truUy detected in the Ks-ba,nd image 
of COSMOS (McCracken et al.), although their associ- 
ated signal to noise is very low. As previously reported 



some of them also have counterparts at 3.6 /im in the 
IRAC data. These "by-eye" detections gave us further 
confidence in the reliability of their identification in the 
i+-band catalog of 109. 

To summarize, we derived photometric redshifts for 
~ 92% objects of the initial 24 jim catalog selected above 
60 /iJy. Sources with counterparts in the Kg-hand but no 
identification at i^g < 26.5 mag account for a fraction of 
5%, while 3% of the catalog correspond to 24 /im sources 
with no identification in the optical nor in the near-IR. 
If we rather consider a selection above 80/iJy, we find 
a total of 29 410 sources in the initial MIPS-24 /im cat- 
alog and 94% of these objects are identified with a red- 
shift. Independently of their possible association with a 
i^g-band or 3.6 /xm counterpart, the 24 /xm sources with 
no redshift are very faint at optical wavelengths. Pre- 
vious studies have shown that they mostly correspond 
to highly-obscured luminous galaxies at 1 < z < 3 rather 
than low redshift sources populating the faint end of the 
luminosity function (e.g., Houck et al. 2005; Dey et al. 
2008). 

2.3.3. Characterization of the photometric redshift 
uncertainties 

By construction an infrared-selected sample of sources 
is biased toward dusty galaxies. The determination of 
photometric redshifts for these objects is subject to more 
uncertainties if the effect of dust extinction is not prop- 
erly taken into account. Furthermore, a large fraction 
of mid-IR sources in the distant Universe arc associated 
with very faint optical counterparts also characterized 
by less accurate redshifts. The uncertainties affecting 
the photometric redshifts of the COSMOS MIPS-24 /im 
sources should thus be carefully assessed before analyz- 
ing our sample further. 

Ilbert et al. (2009) presented a detailed comparison be- 
tween their photometric redshifts and the spectroscopic 
redshifts of MIPS-24 /im sources determined by Lilly 
et al. (2007) and Kartaltepe et al. (in prep.) with spec- 
troscopic follow-ups of the COSMOS field. Up to 2 ^ 1.5 
109 found median dispersions (Taz/(i+z) as low as 0.01 
and 0.05 for 24 /im sources with 22.5 mag < i^g < 24 mag 
and 24 mag < i^g < 25 mag respectively. For mid-IR 
galaxies at higher redshifts and/or with fainter optical 
counterparts there is however no spectroscopic sample 
currently available in COSMOS to allow such compari- 
son. To characterize the uncertainties of their photomet- 
ric redshifts we relied on their associated Probability Dis- 
tribution Functions PDF(z) as provided by 109. For their 
optically-selected sample Ilbert et al. (2009) showed that 
the 1(7 uncertainties'^^ derived from these PDF(z) pro- 
vide a reliable estimate of the dispersion crAz/(i+z) ob- 
tained from the comparison with spectroscopic redshifts. 
In Figure 2 we show these la uncertainties restricted 
to the MIPS-24 /im sources as a function of their pho- 
tometric redshifts, along with their median value calcu- 
lated in three different bins of i+-band magnitude. Up to 

Here we define the positive 1(t+ and negative Icr^ photo- 
metric redshift uneertainties using the relation X'^C^bcst + lo"^) = 
X'^(zbcst ^ lo"^ ) = min(x^(2)) + l, with Zbcst the most Hkely photo- 
metric redshift estimated from the minimum of the merit function 
X^(z). For each object the final la uncertainty is taken as the 
maximum between 1it+ and la~ . 
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the highest redshift bins characterizing our sample (i.e., 
z < 3) we see that the photometric redshifts are still quite 
robust up to < 25 mag (lcr<0.15). At fainter mag- 
nitudes, the uncertainties are somehow larger but the 
dispersions are still as low as la < 0.25 up to z ~ 2 where 
most of the faint sources in our sample are located (see 
Sect. 3). 



0.1 



■ i* < S4 (AB) 

■ 24 < i* < 25 (AB) 

■ 25 < i* < 26.5 (AB) 




1.0 1.5 

Redshift 



Fig. 2. — Photometric redshift uncertainties (~ 1 cr) of the 
COSMOS 24 /xm-selected sources with S24^m > 80 ^Jy (dots), es- 
timated using the probability distribution functions of Ilbert et al. 
(2009) and represented as a function of photometric redshift. The 
sohd hnes represent the median values calculated in three bins 
of i!+-band magnitude (open diamonds: i+<24mag; asterisks: 
24 mag <!"'"< 25 mag; open squares: 25 mag <«"'"< 26.5 mag) and 
using at least 50 sources per bin of redshift. The associated error 
bars were computed from the median absolute deviation of each 
underlying distribution. [See the electronic edition of the Journal 
for a color version of this figure.] 
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Fig. 3. — Fraction of 24 /xm sources with S24(im > 80 /xjy and 
characterized by a redshift probability distribution function with 
two different peaks, illustrated as a function of redshift for the 
same bins of i+-band magnitude as those considered in Figure 2 
(top panel). Each value was estimated using at least a total of 
10 sources per bin of redshift and i+-band magnitude. Error bars 
were computed from Poissonian uncertainties. The bottom panel 
shows the redshift distributions of the 24 ^m sources associated 
with counterparts at 24 mag < j+< 25 mag (light grey shaded re- 
gion) and 25 mag < j+< 26.5 mag (dark grey shaded region). Note 
the very small fraction of outliers up to ^ ~ 2.5. The increasing 
trend observed at i+ > 24 mag and 0.8 ^ z ^ 1.3 affects only a rel- 
atively small number of sources. [See the electronic edition of the 
Journal for a color version of this figure,] 



The larger uncertainties that may affect the photomet- 
ric redshifts of faint optical sources can be produced by 
the presence of multiple peaks in their PDF(z). To quan- 
tify the contribution of these potential catastrophic fail- 
ures we considered the fraction of MIPS-24/zm sources 
with a PDF(z) showing a secondary local maximum 
stronger than 10% of the value of its first peak. These 
fractions are illustrated as a function of redshift and for 
different bins of optical magnitudes in Figure 3. They 
show that up to z ~ 2.5 the number of possible outliers 
is negligible (< 5%) for most of our sample. Fractions 
reaching ~ 25% may be observed at i^g > 24 mag and 
0-8 < 1.3, yet these estimates only affect a very small 
number of sources (see bottom panel of Figure 3). 



3. DIFFERENTIAL SOURCE NUMBER COUNTS 

The differential 2Afim source number counts normal- 
ized to the euclidean slope and derived from our COS- 
MOS data are shown in Figure 4. They are compared 
with the counts obtained by Papovich et al. (2004) based 
on the MIPS Guaranteed Time Observer programs and 
by Chary et al. (2004) using deep observations in the 
ELAIS-Nl field. Within the uncertainties our results 
show a quite good agreement with these two other sur- 
veys, and the variations between the three determina- 



tions may result from the effect of cosmic variance be- 
tween the different fields observed with Spitzer. At fiuxes 
brighter than ~ 3 mJy though, the catalog that we de- 
rived outside of the optically-masked regions of COS- 
MOS underestimates the density of sources compared to 
the counts derived over the whole area of the field. This 
is not surprising given that many of the 24 /zm bright 
objects are associated with optically-bright counterparts 
in these masked regions. Such objects are either stars 
or very low redshift galaxies, which should not affect 
too much our analysis of the general extragalactic 24 /im 
source population. 

Similar to what was observed in other Spitzer 24 fim 
deep surveys (Marleau et al. 2004; Papovich et al. 2004; 
Chary et al. 2004), the counts measured from our COS- 
MOS data show a prominent break at flux densities 
S24:fim 0.2-0.3 mJy. Our source extraction is highly re- 
liable at these flux levels and this feature can not be 
artificially produced by an incompleteness of our 24 /im 
catalog. Such a rapid "turn-over" in the cuclidean- 
normalized source number counts implies their rapid con- 
vergence at fainter fiuxes (<S'24pm ^ 0.1 mJy). Hence, a 
large fraction of the 24 fiin background must be pro- 
duced by sources readily detected in our COSMOS MIPS 
data. For instance Papovich et al. (2004) assumed a 
straight extrapolation of the faint-end slope of the counts 
measured down to 824^1x1 ~ 60 /iJy and they argued that 
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Fig. 4. — DifFerential 24 source number counts derived from 
the COSMOS field outside of the optically-masked regions ('plus' 
symbols). The grey shaded area illustrates the effect of adding a 
systematic uncertainty of it 5% to the absolute 24 /im flux calibra- 
tion. The counts derived for the full coverage of COSMOS (open 
diamonds) as well as the counts from Papovich et al. (2004, aster- 
isks) and Chary et al. (2004, open squares) are shown for compar- 
ison. Our error bars only take into account the statistical Poisson 
uncertainties. Our selection outside the masked regions is mostly 
complete for sources at 80 fi3y < S24fim ^ 3 mjy. The dotted line 
shows our extrapolation of the COSMOS source number counts at 
fluxes below 80/.tJy (see Sect. 5). [See the electronic edition of the 
Journal for a color version of this figure.] 



MIPS-24^m sources detected above 60 ^Jy resolve as 
much as ~ 70% of the 24 /im background intensity. In 
Sect. 5 we will use our photometric rcdshift identifica- 
tions in COSMOS to explore the build-up of this mid-IR 
background light as a function of cosmic time. 

Finally, we note that the persistent "shoulder" appear- 
ing at S2ip,m ~ 3 mJy in the normalized counts that were 
derived from large area surveys such as the Spitzer Wide- 
Area Infrared Extragalactic Survey (SWIRE, Shupe et al. 
2008) is also visible in COSMOS (although our data alone 
can not formally exclude an effect of cosmic variance) . To 
our knowledge this characteristics of the 24 /im counts 
has not been explained so far. 



4. OBSERVED REDSHIFT DISTRIBUTIONS 

We now explore the redshift distribution of the COS- 
MOS MIPS-24/im sources based on the identifications 
of their photometric redshifts described in Sect. 2.3. We 
define H24{z,i^^, Sf-'^'^) the redshift distribution of the 
24 /.tm-selected sources with S24:p,m > Sf^"^ and associ- 
ated with optical counterparts at i+ < ijt^^, and we note 
'^opt(-z, iJim) redshift distribution of all the COS- 
MOS optically-selected sources with i+ < ijt^. Given 
the relatively small number of AGN-dominated sources 
in the general population of optically-selected galaxies^^, 
Ti-optiz, *ii,„) was computed directly using the photomet- 
ric redshift catalog of Ilbert et al. (2009). In Figure 5 

^® Only 0.8% of optical sources with i+ < 24 mag are detected 
in the XMM observations of COSMOS down to Fq 5_2keV ~ 5 X 
10^^^ ergcm^^s"^. 



we show the redshift distributions 7^24(2, ijtjjj, Sf^™) cal- 
culated down to flux limits of 80 ftiy and 300 /iJy at 
24 /xm and for i+-band magnitudes limits of 24, 25 and 
26.5 mag. Sources brighter than i^g = 20mag were ex- 
cluded from all these estimates to minimize the effect of 
large scale structures at low redshift. Our redshift distri- 
bution at < 0.4 should actually be treated with caution 
given the limited comoving volume sampled by our sur- 
vey at these small distances. 

We see that to the depth of our MIPS observa- 
tions (5*24^111 > 80 ^Jy) almost all 24 /xm sources at z < 1 
are associated with optical counterparts brighter than 
i^g~24mag. Our photometric redshifts are very ac- 
curate at these optical magnitudes. Given the statis- 
tics of our sample beyond z^O.5 our redshift distribu- 
tions at 0.5 <z<l should thus be highly robust indepen- 
dently of the 24 /im flux limit considered above 80 /i Jy. 
At z>\, the contribution of fainter sources becomes 
more significant. For instance, 45% and 33% of the 
MIPS-24 /im sources with respectively 5*24^1111 > 80 /iJy 
and 5*24^111 > 300 /i Jy at z > 1 are associated with coun- 
terparts fainter than i^g=24mag. As we showed in 
Sect. 2.3 their photometric redshifts are still quite reli- 
able though (see Figures 2 & 3). 

Our redshift distribution estimates are mainly affected 
by the errors on the COSMOS photometric redshifts, the 
cosmic variance and the errors on the 24 /im flux densi- 
ties. To quantify the impact of all these uncertainties we 
first divided the COSMOS field into 4 different regions 
of similar areas (~0.4deg^) and we derived the possible 
range of redshift distributions in each region with a set 
of 5 000 Monte Carlo simulations. For each of these sim- 
ulations we considered our full 24 /im catalog down to 
'5'24/jm = 60 /iJy and wc modified the 24 /im flux and red- 
shift of each single source according to their 24 /xm flux 
uncertainty and their la photometric redshift dispersion 
characterized in Sect. 2. Assuming a given 24 /zm flux 
limit, an average and independent redshift distribution 
was estimated for the 4 different regions and the disper- 
sion between these 4 estimates was taken as the final un- 
certainty affecting the global 24 //m source redshift distri- 
bution. Our different estimates are reported in Table 1. 
They are also shown in the bottom panels of Figure 5 
for 24 /im flux limits of 80 /i Jy and 300 /i Jy. They sug- 
gest that our global redshift distributions inferred from 
the 1.68 deg^ area observed with MIPS should be robust 
within ~20% up to z~3. 

It is important to note that the redshift distribution 
characterizing the galaxy population selected in the mid- 
IR differs quite significantly from the one observed at op- 
tical wavelengths. Figure 6 illustrates how the optically 
and mid-IR selected sources differ in their occupation of 
the redshift/i+-band magnitude parameter space. For 
instance we find that only ^ 27% of the optically-selected 
galaxies with i < 24 mag are located at z > 1 (Ilbert 
et al. 2009) while this fraction rises to more than 50% 
for the 24 /im sources with S24^m > 80 /i Jy. This large 
proportion of high-redshift mid-IR galaxies is mostly lo- 
cated at 1 < z < 2.5 and it has already been seen in other 
MIPS high redshift surveys (e.g., Houck et al. 2005; Yan 
et al. 2005; Caputi et al. 2006; Papovich et al. 2007; De- 
sai et al. 2008). It results primarily from the strong evo- 
lution that infrared-luminous galaxies have undergone 
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Fig. 5. 



Redshift distributions T-i24{z,i'^^, Sf^'^) of 24 sources with Sf^'^ = 80fi3y (left) and 5^^™=300^iJy (right). The top 

panels represent the redshift histograms determined for different cuts of i Jg— band magnitude as indicated by each legend. The contribution 

of objects with no redshift (i^g > 26.5 mag) is represented as a single bin on the right hand side of the diagrams. The bottom panels show 

our final and best estimates of 7^24 (-^i •S'fj^™) along with their corresponding la uncertainties (white solid lines). They were derived from 
the range of possible redshift distributions defined by the dotted lines and obtained after accounting for the effect of cosmic variance 
and Poissonian uncertainties. The variance produced by large scale structures was quantified by computing the redshift distributions 
in 4 different '^0.4 deg^ regions of COSMOS and with Monte-Carlo simulations accounting for the flux and redshift uncertainties of the 



MIPS sources (grey shaded areas). Galaxies brighter than 



: 20 mag were excluded from all samples to minimize the cosmic variance 



fluctuations at low redshift. [See the electronic edition of the Journal for a color version of this figure. 



with time (e.g., Blain et al. 1999; Franceschini et al. 
2001; Chary & Elbaz 2001). Given the broad-band aro- 
matic features commonly observed at 6 /ini < A < 20 /im 
in the spectra of star-forming galaxies it may also be due 
to galaxy fc-corrections boosting the mid-IR detectabil- 
ity of sources at 1.5<z<2.5 when these features are 
redshifted into the MIPS-24/im bandpass (e.g., Lagache 
et al. 2004). 

To better characterize this large contribution of mid-IR 
selected sources at high redshift, we considered the rela- 
tive distribution ^24(2,1^^, Sf^"^)/Hopt{z,i^i„^)- First, 
this quantity allows us to check if any feature of the 
MIPS-24 /im galaxy redshift distribution does not arti- 
ficially result from a bias already present in the under- 
lying optically-selected catalog of photometric redshifts. 
Besides, it can be seen as the fraction of optical sources 
with z+ < and with a mid-IR detection above a 24 /im 

fiux limit •S'j^^™. Our different estimates of 7^24 /'^opt 
are illustrated in the bottom panels of Figure 6 for the 
same 24 /.tm flux and i~'"-band magnitude limits as the 
ones considered previously. We see that independently 
of these limits the fraction of optical sources detected in 
the mid-IR progressively rises with redshift, which must 
be related to the well-known increase of star formation 
density with lookback time. The fraction peaks at z 2 
before gradually declining at higher redshifts, probably 
because of the decreasing number of 24 /im detections at 



More interestingly though, we observe a substantial 
drop of H24/Wopt between zr^l and z 1.3 followed by 
a steep increase from z~1.3 and z^2. This effect is 
clearly apparent even with a cut at i+<24mag. There- 
fore it can not originate from systematics affecting our 
photometric redshifts. It is rather due to the break char- 
acterizing the 24 /im source redshift distributions around 
z ~ 1.3 and that does not appear in the smooth distribu- 
tions observed at optical wavelengths (e.g., 109). Note 
that this break also seems to be more prominent as the 
24 /im flux limit increases. The redshift distribution de- 
rived at S24fim > 300 /i Jy looks even like a two-component 
function with a main and narrow peak of sources at z ~ 1 
followed by a secondary but broader peak at 1.5 <z< 2.5 
(e.g.. Figure 5, right panel). This characteristic of the 
24 /im redshift distribution could be artificially produced 
by a relative increase of mid-IR detections at z ~ 2 due 
to the rest-frame 7.7 /im aromatic feature enhancing the 
apparent 24 /im flux of star- forming galaxies when it is 
redshifted into the MIPS-24 /im bandpass. On the other 
hand it may also reflect a smaller rate of detections at 
z ^ 1.3 given the possible effect of the 9.7 fim silicate ab- 
sorption. If this break is truUy produced by fc-correction 
effects, the stronger effect observed at brighter 24 /im flux 
might also suggest more prominent features in the more 
luminous mid-IR selected galaxies. 
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TABLE 1 

Redshift distributions of the cosmos MIPS-24 fj,u selected sources 



Redshift 


N(z, S24>0.08mJy) 


N(z, S24>0.15mJy) 


N(z, S24 > 0.3 mj 


0.0<2<0.2 


679 ±91 


310 ±63 


110 ±20 


0.2<2<0.4 


4722 ± 203 


2356 ± 64 


697 ±62 


0.4<2<0.6 


6317 ±275 


3296 ± 71 


997 ±44 


0.6<2<0.8 


9991 ± 1328 


4964 ± 796 


1604 ± 353 


0.8<2<1.0 


14055 ± 2001 


7144 ±1193 


2148 ±377 


1.0<2< 1.2 


8981 ± 904 


3743 ± 471 


865 ±117 


1.2<z<lA 


5158 ±465 


1704 ± 167 


367 ±54 


1.4<z< 1.6 


6360 ± 144 


2356 ± 155 


531 ± 78 


1.6<z< 1.8 


4672 ± 288 


1993 ±117 


452 ±67 


1.8<2<2.0 


4735 ± 264 


2167 ±203 


497 ±30 


2.0<2<2.2 


3753 ± 203 


1683 ± 150 


411 ±62 


2.2<2<2.4 


1933 ±65 


815 ±52 


205 ± 13 


2.4<2<2.6 


1423 ±93 


598 ± 51 


167 ±27 


2.6<2<2.8 


978 ±110 


406 ± 54 


126 ± 31 


2.8<2<3.0 


535 ± 89 


222 ± 32 


84 ± 21 


3.0<z<3.2 


335 ± 34 


125 ±8 


49 ±6 


3.2<2<3.4 


203 ± 30 


103 ± 17 


43 ±5 


3.4<2<3.6 


138 ±7 


67 ±14 


36 ±18 


3.6<2<3.8 


73 ±7 


42 ±8 


18 ±7 


3.8<2<4.0 


53 ±8 


31±9 


16 ±5 


4.0 < z < 1.2 


46 ± 10 


34 ± 7 


16 ± 3 



(^) The distributions N(z) refer to the number of sources per redshift unit and per square 
degree. The la error bars account for cosmic variance, Poisson noise, flux and redshift 
uncertainties as described in Sect. 4. 



5. THE BUILD-UP OF THE MID-INFRARED 
BACKGROUND ACROSS COSMIC AGES 

The Cosmic Infrared Background (CIB) refers to the 
extragalactic background hght (EEL) detected between 
~ 8 fim and ~ 1000 /itm. It can be seen as a fossil record 
of all the radiations emitted at IR wavelengths along 
the process of galaxy formation (e.g., Elbaz ct al. 2002). 
Its integrated energy power is comparable to that of the 
background produced by the stellar component of galax- 
ies in the optical/near-IR wavelength range (Dole et al. 
2006). Understanding how the different populations of 
distant sources contributed to the build-up of the CIB 
as a function of redshift can thus provide significant 
constraints to our general understanding of the cosmic 
growth of structures. 

Using stacking analysis at 70/Ltm and 160 /xm Dole 
et al. (2006) showed that 24 /im--sclcctcd sources with 
'S'24pm > 60 /iJy contribute more than 70% of the CIB. 
The large statistics available in the COSMOS field and 
the high rate (~95%) of redshift identifications in our 
24 nm catalog give us therefore a unique opportunity 
to resolve the different contributions of galaxies to the 
build-up of a large fraction of the EBL. In this current 
paper we only explore the history of the monochromatic 
background light at the wavelength of our MIPS observa- 
tions, using the differential 24 /um source number counts 
as a function of lookback time. Resolving the history of 
mid-IR source number counts already provides tight con- 
straints on the different models of galaxy formation, espe- 
cially when they correctly reproduce integrated measures 
like the total counts or the diffuse background emission 
but differ in their predictions of the flux/redshift distri- 
bution of galaxies responsible for these quantities. This 
issue will actually be discussed in Sect. 6 where wc will 
compare our data with the predictions of different mod- 



els to quantify their reliability in the description of IR 
galaxy evolution. We defer to a subsequent publication 
the analysis of the total background intensity produced 
by 24 /Ltm-selected galaxies over the entire spectral range 
of the CIB. 

To characterize this decomposition of the 24 /im 
counts as a function of redshift we first defined 

dN{zii-cn,S24ij.m)/dS24:t^m and dN{i^^,S24f_im)/dS24^m 

the differential counts of the 24 /xm sources identified up 
to a given redshift znm or up to a given z+-band mag- 
nitude i'^^, respectively. In Figure 7 we illustrate these 
quantities computed for ziim = 0.5, 1, 1.5, 2 and 2.5 (left 
panel) and for i+^ = 21, 22, 23, 24, 25 and 26.5 mag 
(right panel). For comparison we also show the total 
number counts derived from the COSMOS field and dis- 
cussed in Sect. 3 (sec also Figure 4). We sec that the 
brightest 24 /xm sources are mostly associated with low 
redshift galaxies characterized by bright optical counter- 
parts. However a non-negligible fraction of the MIPS de- 
tections at 524pm ^ 1 mJy are also located at much larger 
distances. Albeit rare in space density, these sources 
compose the bright end of the IR luminosity function at 
high redshift and an accurate estimate of their contribu- 
tion can only be determined with large area surveys. At 
fainter fluxes, we also note that a substantial fraction of 
the 24 /im sources are associated with faint and distant 
optical counterparts at ij3>24mag. This conclusion 
reinforces that discussed above for the redshift distribu- 
tions. It emphasizes once again the difficulties to achieve 
a complete characterization of the 24 /xm-selected galaxy 
population with current optical/near-IR follow-ups, and 
in particular to obtain a full census of the sources at 
5*24/^01 ^ 0.5 mjy making the bulk of the mid-IR back- 
ground light. At such faint optical magnitudes, spec- 
troscopy has been very challenging so far while photomet- 
ric redshifts also need to be carefully treated given the 
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Fig. 6. — The top panel shows the distribution of i"*"- 
band magnitudes as a function of redshift computed for the 
COSMOS j+-band selected galaxies (Ilbert et al. 2009, 
shaded region, darkness increasing with source density) and 
the 24 /im sources with 524(im > 80 /xjy (white contours). The 
two other panels illustrate the normalized relative distribu- 

•+ (rr24^m\ /rt/ / 

*lim''^lim )/ nopt{Z,l^.^^ .jj^- 



tions W24(^, 
diamonds). 



lim ' 



: 25 mag (asterisks) and jjt^ : 



: 24 mag (open 
= 26.5 mag (open 



squares), with Wopt(2, ^J^^) the redshift distribution of the 
optically-selected sources with < ijt"^ derived by Ilbert et al. 
(2009). They were computed for Sf;^™ of 0.3 mjy and 0.08 mjy 
(middle and bottom panels, respectively). The fraction of optical 
sources detected at 24 ^m is clearly rising up to z 2, although a 
prominent break is visible at 2:~ 1.3. [See the electronic edition of 
the Journal for a color version of this figure.] 



rapid increase of their uncertainties beyond I ^25 mag 
(e.g., 109). 

To also account for the contribution of galaxies fainter 
than our 24 /zm detection limit we extrapolated the 24 fim 
number counts at faint fluxes following the same ap- 
proach as considered by Papovich et al. (2004). We used 
power-law functions dN/dS24:p,raO^S2ip,ui~°' for the full 
24 /im catalog and for the different sub-samples of 24 /im 
sources selected per bin of redshift or i^-band magni- 
tude as described above. For each of these extrapola- 
tions the power-law index a was only constrained with 
the counts observed above S24p,ra = 80 fiJy to minimize 
the impact of the incompleteness of our source extrac- 
tion below this limit. For the total counts we derived 
a~1.5. This is similar to the results obtained by Pa- 
povich et al. (2004), Chary et al. (2004) and Rodighiero 
et al. (2006) who measured indexes between 1.5 and 1.6 
using sources detected down to ~ 30 /i Jy in the ELAIS- 
Nl field. This good agreement suggests that the extrap- 
olations obtained for the other sub-samples should also 
provide a fair representation of the flux and redshift dis- 
tributions characterizing the MIPS sources below our de- 
tection limit at 24 /im. 

Finally, we integrated the total extrapolated number 



counts to estimate the total 24 /im background light (re- 
ferred as Btot hereafter). We also integrated the extrap- 
olated functions dN/ (iS'24/jm x S24fim for each sub-sample 
considered above to define 



/>oo 

^24/xm(^lim, S24i.im) = / dN{ziu^, S) / dS X S X dS 

(1) 

and 



S24Mm(«lt„, ^24^™) = / rfA^(*L, S)/dS xSxdS 

(2) 

as the background light produced above a 24 /im flux 
S24fj.m and up to the redshift (Eq. 1) or by sources 
brighter than the i+-band magnitude i^^^^ (Eq. 2)^^. 

These two quantities are respectively illustrated in 
Figures 8 & 9 for 24 /im fluxes of 0.06, 0.08, 0.15, 0.3 and 
1 mJy as well as for S2A^m = 0. They were normalized to 
the total 24 /zm background intensity _Btot determined 
from the extrapolation of the total source number 
counts. They can thus be read as the fraction of the 
24 /im background resolved by MlPS-selected sources as 
a function of their flux 6*24^™, their redshift and their 
optical magnitude. The bottom panel of Figure 9 also 
illustrates the differential 24 /zm background intensity 
produced as a function of redshift and computed for 
the same 24 /tm fluxes as in the upper panel. Note 
that because of the magnitude limit imposed by our 
optical counterpart identiflcations i?24;im(ii|,„, 'S'24;^m) 
was only computed up to i+ = 26.5 mag. Further- 
more our determination of i324Aim(-ziim, <S'24pm) is 
necessarily underestimated at high redshift given the 
incompleteness of our identifications at faint fiuxes. 
However we believe that these two effects should af- 
fect our results by at most 10% given our estimate of 



B 



24/jm 



(«lim' S2 



'lim' '^24^tm)/-Btot ~ ^24/jm (■^lii 



for 5. 



24;im 



= 0, 



'lim 



— 26.5 and zn, 



, S24iim)/ Btot ' 

00. 



'0.9 



While mid-IR shallow surveys (e.g., 524^™ > 0.3 mJy) 
only resolve up to ~30% of i3tot the analysis of Fig- 
ures 8 & 9 reveals that the MIPS-24/im sources brighter 
than S24^im ^ 60 /i Jy resolve up to ^ 70% of the total 
24 /im background. This estimate is consistent with that 
of Papovich et al. (2004) but it is slightly higher than 
the results obtained by Rodighiero et al. (2006) in the 
ELAIS-Nl field. Since a non-negligible fraction of the 
background is produced by bright sources (e.g., > 10% 
at S24fj.m ^ 1 mJy) this disagreement could be simply ex- 
plained by the small area of the sky covered in their 
analysis (~ 185 arcmin^) and the small number of bright 
24 /im detections in their data. Furthermore up to ~ 10% 
of the 24 /tm background originates from optical sources 
fainter than the limit of our photometric redshift cata- 
log (i+ = 26.5 mag), but more than 50% of Btot is due to 
sources detected above 80 /iJy at z < 1.5 and associated 
with optical counterparts brighter than i+ = 24 mag. Al- 
though their contribution to the integrated CIB would 

The spectral slope of the counts at faint fluxes guarantees the 
convergence of the integrals when S2Aiim — > whatever the choice 
of the redshift or the i+-band magnitude limit. 
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Fig. 7. — Differential MIPS-24 /im number counts along with the contributions of 24 sources computed for different cuts in redshift 
(left panel) or j+-band magnitude (right panel) and as shown by the shaded regions. The thick black curve represents the total number 
counts derived from the whole 24 /xm catalog. The bin denoted as "All z" in the left panel corresponds to all 24 /^m sources identified with 
a redshift in our survey. The corresponding fractions of sources observed in the different bins of selection are illustrated as a function of 
24 /im flux density in the bottom panels. Note the increasing contribution of optically-faint (ijg > 24 mag) and high-redshift galaxies to 
the source number counts at faint 24 /xm fluxes. [See the electronic edition of the Journal for a color version of this figure.] 



require a careful analysis of their SED in the far-IR, 
our picture is therefore broadly consistent with previ- 
ous claims supporting a large contribution of infrared- 
luminous galaxies at0.5<2:<1.5to the EBL (e.g., Elbaz 
et al. 2002; Dole et al. 2006). Conversely, we also find 
that up to ~ 30% of the total 24 /im background must 
be produced by optically-faint (i+ > 24 mag) galaxies at 
z>l. The exact redshift distribution of these sources 
can not be constrained with our data given our 24 /im 
flux limit of 80 iiJy as well as the incompleteness of our 
redshift determinations at faint 24 /xm fluxes. However, 
the flattening of i?24/jm(-ziim, 5'24pm) observed in Figure 9 
suggests that they should be mostly located at ^^z<3. 

Finally, the bottom panel of Figure 9 reveals that the 
relative contribution of galaxies to the 24 /im background 
increases up to ~ 1 and rapidly declines at higher red- 
shift. Although this trend clearly depends on our ex- 
trapolations of the counts below 80 /iJy it is again con- 
sistent with the picture suggested by phenomenological 
models of IR galaxy evolution (e.g. Chary & Elbaz 2001; 
Franceschini et al. 2001) and by the analysis of the his- 
tory of the CIB recently performed with other methods 
(e.g.. Dole et al. 2006). Nonetheless, the contribution of 
galaxies with 6*24^111 > 80 /iJy seems to be roughly con- 
stant up to 2:~1. Hence, the "peak" observed at z ~ 1 
must be produced by a strong evolving (albeit not dom- 
inant) contribution of galaxies fainter than 80 /iJy and 
with rather low infrared luminosities. 



6. COMPARISON WITH MODEL PREDICTIONS 
6.1. Semi- analytical models 
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Fig. 8. — Fraction of the total 24 ^m background produced 
by sources brighter than the j+-band magnitude (black thick 
solid line, assuming an extrapolation of the counts at the faint end), 
along with the fractions resolved above 24 /im fluxes of 0.06, 0.08, 
0.15, 0.3 and 1 mJy (solid lines, increasing thickness corresponding 
to fainter flux limits). For comparison the blue-dashed (resp. red 
dash-dotted) line shows how the i+-band (resp. Xs-band) sources 
contribute to the fraction of the 24 /im background produced above 
S'24/jm = 80 /ijy and directly measured in our data. 



Constraining the physics of galaxy formation with 
semi-analytical models has become over the last decade 
one of the major challenges of modern cosmology. Based 
on the now-popular cold dark matter (CDM) paradigm 
of structure formation (Springel et al. 2005) these models 
compute the evolution of the stellar content of galaxies 
from the hierarchical assembly of dark matter haloes, us- 
ing detailed simulations to account for a variety of physi- 
cal processes such as the radiative cooling and the shock- 
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Fig. 9. — Top: Fraction of the total 24 fim background pro- 
duced by the COSMOS-selected i+ <26.5 mag sources below the 
redshift z^^^ (black thick solid line, assuming an extrapolation of 
the counts at the faint end) , along with the fractions resolved above 
2Afim fluxes of 0.06, 0.08, 0.15, 0.3 and 1 mjy. For comparison 
the blue-dashed line shows how the identified 24 fim sources con- 
tribute to the fraction of the 24 ^m background produced above 
'S'24;jm = 80 fi3y and directly measured in our data. Uncertain- 
ties due to the incompleteness of our redshift identification beyond 
2 ~ 1 should be at most 10%. Bottom: Differential background 
intensity (arbitrary units) produced as a function of redshift and 
computed with the same selection criteria as in the top panel. 



heating of gas, the energy feedback from supernovae and 
active nuclei, the chemical enrichment of the interstel- 
lar medium and galaxy mergers (e.g., Cole et al. 2000; 
Bower et al. 2006; Croton et al. 2006; Cattaneo et al. 
2006). Although many quantitative aspects of galaxy 
evolution are still not correctly predicted with these sce- 
narios they have recently been successful in reproduc- 
ing several key properties of present-day galaxies such 
as their large-scale clustering or the bimodality of their 
optical color distribution (Blaizot et al. 2006; Somerville 
et al. 2008). 

Because of the additional difficulty to properly account 
for dust extinction, the determination of accurate predic- 
tions of the deep IR Universe with these models has been 
however barely explored so far. Simulating the evolution 
of galaxy IR properties requires not only a correct es- 
timate of the history of the physical processes powering 
their energy output over relatively short time scales (e.g., 
star- forming activity, nuclear accretion), but an accurate 
treatment of their reprocessing of UV/optical radiation 
at longer wavelengths is also critical. Since the launch 
of Spitzer Lacey et al. (2008) published the only single 
study (to our knowledge) describing some of the expected 
properties of 24 /im-selected galaxies with the use of a 
semi-analytical model based on the hierarchical growth of 
dark matter haloes^^. Their simulations were performed 
with the GALFORM galaxy formation model (Cole et al. 
2000). They used the GRASIL spectro-photometric code 
(Silva et al. 1998) to calculate SEDs and to account for 
the reprocessing of radiation by dust, which allows for 

Silva et al. (2005) derived similar predictions based on the 
semi-analytical model of Granato et al. (2004), but this model does 
not include galaxy merging. 



predictions of the multi-wavelength properties of all sim- 
ulated galaxies as a function of cosmic time. 

Figure 10 shows their model predictions for the evolu- 
tion of the galaxy median redshift as a function of 24 fim 
flux, along with their estimates of the 10*'' and 90**^ per- 
centiles characterizing the associated redshift distribu- 
tions. They are compared with the same quantities esti- 
mated from our COSMOS data, where the reported error 
bars account for the combination of flux and redshift un- 
certainties, the cosmic variance and the Poisson noise as 
described in Sect. 4. Interestingly enough, we note that 
the model of Lacey et al. (2008) shows a statistically- 
significant overestimate of the median redshift of galaxies 
at all fluxes probed by our survey. This overestimate is 
particularly noticeable at fluxes below ^0.5mJy. Given 
that their predictions agree reasonably well with the to- 
tal source number counts observed at such flux levels (see 
their figure 2), it implies an excess of faint sources pre- 
dicted above z ~ 1 as well as a lack of sources predicted at 
lower redshifts. This disagreement was already noticed 
by Lacey et al. (2008) from the comparison between their 
model and luminosity functions previously derived from 
Spitzer data. It is also clearly apparent in Figure 11, 
where we compare the COSMOS redshift distribution 
measured above 824^,171 — 83 /iJy with the result of the 
simulations that they obtained for the same 24 fim flux 
limit. On the other hand we also note that the 90**^ per- 
centile of the distribution observed in COSMOS at bright 
fluxes seems to extend up to higher redshifts compared 
to the model predictions. This could be interpreted as 
an underestimate of the bright end of the simulated IR 
luminosity function at 1.5 <z< 2.5. It could also be due 
to an excess of sources at much lower redshifts since the 
model of Lacey et al. (2008) overestimates the differential 
24 fiia source number counts above ~ 1 mJy. We will in- 
vestigate these different possible interpretations in future 
papers where we will present the IR luminosity functions 
derived from our MIPS imaging of the COSMOS field. 

In Figures 12 and 13 we now confront the observed and 
the predicted distributions of 2Afim sources as a func- 
tion of their ncar-IR and optical photometry. Following 
the same selection criteria as those considered by Lacey 
et al. (2008) we restrict our sample to sources brighter 
than 5*24^111 — 83 //Jy. We split these sources in different 
bins of Ks~ha,nd and r+-band magnitude and we com- 
pare the redshift distributions observed in these bins with 
the corresponding predictions derived from the simula- 
tions^^. The range of distributions shown for the COS- 
MOS data (grey shaded regions) accounts for our best 
estimate of cosmic variance, fiux and redshift uncertain- 
ties (see Sect. 4). The photometric errors affecting the 
r+-band and Ks-hand data were not taken into account 
but their effect should be limited given the adopted bins 
of magnitudes and the depth of the COSMOS observa- 
tions (r+ < 26.8 mag AB and _R's<23.8mag AB at 5a: 
Capak et al. 2007, McCracken et al. submitted). 

As already noted by Lacey et al. (2008), our fig- 
ures 12 & 13 show that the excess of simulated sources 

Lacey et al. (2008) provide their predicted distributions using 
magnitudes in the Vega system. We thus transformed our AB 

magnitudes assuming '"yoga = ''ab ~ ^-^^ '^'^'^ ^s.Vega = ■'^'s.AB - 
1.84. The transformation between our set of filters (r"*" & Kg) and 
that used by Lacey et al. (R & K) was however neglected. 
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Fig. 10. — Comparison between our COSMOS fiux/redshift dis- 
tributions (thin lines and open symbols) and the predictions from 
the semi-analytical model of Lacey et al. (2008, thick curves, no 
symbol). For each bin of 2Afim flux the solid lines correspond to 
the median redshift while the dashed and dash-dotted lines repre- 
sent the 10*^ and 90**^ percentiles of the associated distribution. 
Our error bars account for the uncertainties due to cosmic variance 
and Poisson noise as described in Sect. 4. The shaded region illus- 
trates the parameter space populated by the MIPS-24 /im sources 
between the 5^^ and the 95^^ percentiles of their redshift distri- 
bution. This comparison reveals a systematic overestimate of the 
predicted median redshift at all 24 fim fluxes. [See the electronic 
edition of the Journal for a color version of this figure.] 
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Fig. 11.— Total redshift distribution of the COSMOS 24 /xm 
sources with S24(jm > 83 /.tjy (grey shaded region), compared with 
the model predictions of Lacey et al. (2008, solid line and filled 
histogram). The observed COSMOS distribution accounts for the 
uncertainties due to cosmic variance and Poisson noise as described 
in Sect. 4. Note the significant excess of sources predicted at z ~ 2 
as well as the lack of 24 /im simulated galaxies at z ^1. 



appearing in their model occurs mostly around z~2. 
It originates from galaxies fainter than /C^^ 21.5 mag 
(AB) while most of the 24 /im sources in our selection 
(<5'24/.im > 83 n3y) are typically brighter than this mag- 
nitude limit (see also, e.g., Papovich et al. 2007). On 
the other hand our comparison in the r+~-band reveals 
that these simulated sources are not particularly fainter 
than the typical 24 fim high redshift sources at opti- 
cal wavelengths. In fact the redshift distribution of the 
24 fim sources with faint optical fluxes seems to be rea- 
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Fig. 12. — Redshift distributions of COSMOS 24 fim sources with 
S'24pm > 83 ^Jy decomposed into 3 bins of i^s-band Vega magni- 
tude (grey shaded regions), compared with the model predictions of 
Lacey et al. (2008, solid lines and filled histograms). The observed 
distributions account for the uncertainties due to cosmic variance 
and Poisson noise as described in Sect. 4. Note the significant ex- 
cess of sources predicted at z ~ 2 and mostly due to galaxies fainter 
than Ks ~ 21.5 mag, as well as the lack of 24 /ira simulated sources 
with bright near-IR counterparts. 



sonably well predicted (although the uncertainties due 
to the incompleteness of our redshift identification be- 
yond z^g = 26.5 mag prevents reaching definite conclu- 
sions). The excess of galaxies in the simulations orig- 
inates from sources lying among the brightest optical 
counterparts associated with the 24 /im galaxy popula- 
tion at z ~ 2 (i.e., r"*" < 25 mag Vega) and it may thus be 
related to a too large number of blue galaxies predicted 
by the models. 

Similarly we find that the lack of simulated galax- 
ies noted by Lacey et al. (2008) at z ~ 1 is related to 
an underestimate of sources with relatively bright near- 
IR counterparts, while the predicted number of galaxies 
with fainter i^g-band luminosities {Kg > 20.5 mag Vega) 
is clearly too large at this redshift. Given the slightly bet- 
ter agreement found at optical wavelengths, these results 
reveal again some discrepancies not only in the absolute 
optical/near-IR magnitudes of the predicted galaxies but 
also in their colors. We emphasize that this strong lack 
of 24 ^m sources with bright near-IR luminosities (i.e.. 
Kg < 20.5 mag Vega) is found at all redshifts from z ~ 0.5 
to z~2.5. Since the ^^^"band emission of these distant 
sources should probe the rest-frame properties of their 
evolved stellar populations this may reflect a systematic 
underestimate of the density of massive galaxies with on- 
going star formation in the distant Universe. Unfortu- 
nately though, we do not have access to the mock cat- 
alogs where the simulated 24 /im sources were extracted 
from. More detailed comparisons between the physical 
properties of the observed and the modeled galaxies will 
have to be performed to better understand the origin of 
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Fig. 13.— Redshift distributions of COSMOS 24 /xm sources 
with S'24(jm > 83 /ijy decomposed into 4 different bins of r+-band 
Vega magnitude (grey shaded regions), compared with the model 
predictions of Lacey et al. (2008, solid lines and filled histograms). 
The observed distributions account for the uncertainties due to 
cosmic variance and Poisson noise as described in Sect. 4. The 
predictions show a statistically better agreement than in the Kg- 
band (Fig. 12), although strong discrepancies can still be noted. 



the discrepancies discussed above. 

6.2. Phenomenological models 

As shown above the understanding of the deep IR Uni- 
verse with semi-analytical simulations still remains an 
important challenge for extragalactic astronomy. In the 
past few years, simpler descriptions accounting for the 
strong evolution of galaxy IR properties and based on 
"phenomenological" models have been proposed in the 
literature. Assuming a backward evolution of the local 
luminosity function and a library of galaxy spectral en- 
ergy distributions, these models are usually constrained 
to reproduce a number of quantities like the Cosmic In- 
frared Background and/or galaxy IR properties such as 
the source number counts and redshift distributions. Al- 
though they can not be used to follow the history of 
the physical ingredients driving the evolution of galaxies 
and the change of their properties with time (e.g., gas 
accretion, stellar mass build-up, interactions, feedback, 
...), they do not require a formalism as complex as the 
one needed in semi-analytical models. They can also be 
useful in planing the overall characteristics (e.g., source 
counts, redshift distribution, ...) of new IR or submil- 
limeter surveys to be performed at yet-unexplored wave- 
lengths. 



All these models can vary from one another in the 
choice of their IR SED library, their decomposition of the 
overall galaxy population in diS'erent sub-groups (e.g., 
quiescent galaxies, starbursts, active nuclei, ...) or the 
parameterization characterizing how the luminosity func- 
tion is allowed to vary with lookback time. In Figure 14 
we present the comparison between our COSMOS 24 /xm 
source redshift distributions and the predictions from 
the phenomenological models of Le Borgne et al. (2009), 
Gruppioni et al. (2005, in prep.) and Lagache et al. 
(2004). Assuming the luminosity-dependent SED library 
from Chary & Elbaz (2001), Le Borgne et al. (2009) have 
used a non-parametric "count inversion" technique with 
no specific assumption on the evolution of the luminos- 
ity function to fit the most-recently updated results from 
Spitzer IR surveys. On the other hand the model pre- 
sented by Gruppioni et al. (2005) did not make use of 
any Spitzer data but was constrained with 15 /zm obser- 
vations carried out with the Infrared Space Observatory . 
It assumes four distinct populations of sources (quiescent 
galaxies, star-forming sources, type 1 and type 2 AGNs), 
each of them evolving independently from the others. 
Initially these populations were characterized by a single 
template SED (Gruppioni et al. 2005), but their scenario 
was recently updated to account for new procedures fit- 
ting individual galaxies with AGN and starburst SEDs 
simultaneously from the optical up to the mid-IR wave- 
lengths (Gruppioni et al., in prep.). Finally, Lagache 
et al. (2004) assume a population of luminous starbursts 
characterized by luminosity-dependent SEDs as well as 
a population of quiescent cold galaxies evolving indepen- 
dently from the star-forming sources. Their model was 
updated from a previous version developed before the 
launch of Spitzer to account for the new constraints im- 
posed at 24 fim by the MIPS cosmological surveys. 

The comparisons between these three models and our 
COSMOS data are illustrated for 24 /im flux limits of 
0.08, 0.15 and 0.3 mJy. The possible range of observed 
redshift distributions shown in the figure was calculated 
as described in Sect. 4 by taking into account the com- 
bined effect of flux and redshift uncertainties as well 
as the Poisson noise and the cosmic variance. We did 
not include any comparison below z '^OA as the small 
comoving volume sampled by our survey at these low 
redshifts would imply large statistical uncertainties. At 
z>l and faint fluxes (e.g., S'24;jm ^ 0.15mJy) we note 
that the model of Le Borgne et al. (2009) shows a rather 
good agreement with our observations, while the models 
from Gruppioni et al. (2005) and Lagache et al. (2004) 
show a significant overestimate of the number of 24 /im 
sources around z^2. As we discussed in Sect. 4 the 
bimodal distribution predicted by these two models is 
mainly produced by /c-correction effects boosting the de- 
tectability of starburst sources when their 7.7 /im aro- 
matic feature is redshifted into the MIPS 24 /im band- 
pass at 1.5<z<2.5. The predicted bimodality is clearly 
not as prominent in the data. This disagreement could 
refiect either the presence of too strong features in the 
SEDs assumed by the models or a too rapid evolution 
of the source space density at the faintest fiuxes typ- 
ically probed by our survey. On the other hand, the 
break observed at 1.3 in the redshift distribution of 
the COSMOS galaxy population at bright fiuxes seems 
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Fig. 14. — Comparison between our COSMOS redshift distributions (grey shaded regions) and the model predictions (solid lines and 
filled histograms) of Le Borgne et al. (2009, left panels), Gruppioni et al. (2005, middle panels) and Lagache et al. (2004, right panels), 
shown for 24 fim flux limits of 0.08, 0.15 and 0.3 mjy. The observed distributions account for the uncertainties due to cosmic variance and 
Poisson noise as described in Sect. 4. The dashed lines in the medium-column panels represent the predictions from a new version of the 
model developed by Gruppioni et al. (in prep., Gruppioni & Pozzi, private communication). 



to be better reproduced by the model of Gruppioni et al. 
than in the predictions from Le Borgne et al. (where the 
break is barely apparent) or Lagache et al. (where the 
break is too prominent due to the strong excess of sources 
expected at z~2). At z^l Gruppioni et al. (2005) pre- 
dict however too many sources at S'24/jin ^ 0.3 mJy while 
the model of Le Borgne et al. may also overestimate the 
luminosity function at 0.4 < z < 0.8. 



7. DISCUSSION 

Because of (i) the unprecedented accuracy of the 
photometric redshifts used in this work, (ii) our very 
high rate of redshift identifications obtained down to 
S24:fim = 80 /iJy, and (iii) the large area covered by our 
survey, our analysis represents the most statistically sig- 
nificant characterization of the redshift distributions of 
the 24 /im source number counts achieved so far. We 
found that the distribution of the 24 /im selected galax- 
ies is characterized by a prominent peak at 2 1 and a 
relatively large fraction of sources at 1 < z < 3, which is 
globally consistent with the distributions of MIPS-24 fim 
sources previously determined with other surveys (Perez- 
Gonzalez et al. 2005; Caputi et al. 2006; Papovich et al. 
2007; Desai et al. 2008). However, specific differences be- 
tween our results and these other studies can be noticed. 
For instance, the redshift distribution obtained by Perez- 



Gonzalez et al. (2005) shows a gradual and rapid decline 
from z 1 up to z ^ 3, while we observe a break at z ~ 1.3 
as well as a larger fraction of sources at z > 1. Similarly, 
our distributions exhibit smoother variations than the 
redshift distribution determined by Caputi et al. (2006) 
in the field of the Chandra Deep Field South (CDFS). 
All these deviations from our results could originate ei- 
ther from less accurate photometric redshifts that were 
used in some of these other studies, or from the effect of 
cosmic variance over small fields of view like the CDFS. 

The large number of 24 /j,m sources at 0.5 < z < 1.5 and 
our analysis of the history of the 24 fj,m background light 
also shows that the GIB must originate predominantly 
from moderately-luminous galaxies at intermediate red- 
shifts rather than highly obscured and luminous dusty 
sources at z>2. This result has already been infered 
from previously-published studies (e.g., Elbaz et al. 2002; 
Dole et al. 2006) but our data provide for the first time 
a large enough statistics to explore with good accuracy 
the build-up of the mid-IR background as a function of 
cosmic time. The characterization of the far-IR proper- 
ties of the 24 /im sources with the already-existing MIPS- 
70/160 /im and MAMB0-1.2mm observations of GOS- 
MOS (Prayer et al., submitted; Bertoldi et al. 2007) or 
with the future coverage of this field with forthcoming fa- 
cilities like Herschel will provide an even more detailed 
picture of the history of the EBL over the whole range 
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of IR wavelengths. 

Yet, in spite of this globally coherent picture that has 
emerged from mid-IR surveys our comparison between 
the COSMOS 24 /xm data and the model predictions of 
Lacey et al. (2008) reveal that the evolution and the 
properties of luminous galaxies in the distant Universe 
are still far from being properly reproduced by the semi- 
analytical simulations of cosmic history. Although the 
rapid evolution of the IR luminosity function up to ^; ^ 2 
is globally well predicted by the model of Lacey et al. 
(2008), discrepancies in the expected source density per 
bin of IR luminosity and in the optical/near-IR proper- 
ties of 24 /nm sources are clearly apparent. 

Interestingly we found that the number of 24 /xm 
sources with bright near-IR counterparts (i.e., 
Kg < 20.5 mag Vega) in the simulations is clearly 
miderestimated at all redshifts beyond It 
implies that the predicted source density of massive 
galaxies with strong on-going activity of star formation 
such as those recently discovered at high redshift with 
Spttzer (e.g., Daddi et al. 2005; Reddy et al. 2006; 
Papovich et al. 2006; Wuyts et al. 2008) is too low. 
Besides we note an overestimate of simulated sources 
with 20.5 mag < iiTs < 21.5 mag (Vega) at z^l, while 
their counterparts in the r+ band seem to be more 
representative of the optical properties of 24 /im sources 
at this redshift. In addition to the overall imdcresti- 
mate of the predicted source density at this redshift, 
we thus conclude that a large fraction of simulated 
galaxies at 2; ~ 1 are too blue with respect to the typical 
24 /xm detected sources. This trend actually recalls the 
discrepancy that we also observed at z ^ 2 where the 
overestimate of predicted sources originates from blue 
galaxies with JCg-band luminosities much fainter than 
the bulk of 24 /xm sources found at this redshift. 

The origin of this excess of blue galaxies with a 24 /xm 
detection above our sensitivity limit is not straight for- 
ward to determine. Assuming that the star formation 
rate (SFR) of individual sources is correctly predicted by 
the simulations, it may be related to an incorrect descrip- 
tion of their spectral energy distributions. As we men- 
tioned in Sect. 4 the broad aromatic features commonly 
observed at 7.7 /im, 11.3 fim and 12.7 /im in the spectra of 
starburst galaxies can enhance their detectability when 
these lines are redshifted into the MIPS-24/um bandpass. 
An overestimate of the strength of the 7.7/xm feature 
could thus result in too many 24 sources predicted at 
z^2. Furthermore the band emission of galaxies at 
such high redshifts relates to their rest-frame UV proper- 
ties, which are obviously very sensitive to dust extinction 
given our selection performed at mid-IR wavelengths. A 
wrong treatment of the obscuration affecting these 24 /im 
sources would thus bias their predicted characteristics at 
optical wavelengths as well as the estimate of their re- 
processed dust emission in the mid-infrared. 

On the other hand, this excess could also originate from 
an overestimate of the galaxy star formation rate itself. 
In the simulations of Lacey et al. (2008) galaxies with 
on-going bursts of star formation account for ^ 50% and 
> 90% the 24 /zm source population at z ~ 1 and z ~ 2 
respectively (the other sources being classified as quies- 
cent galaxies, see their Figure A4). Their mid-IR emis- 
sion directly scales with their star-forming activity, which 



could lead to a too large number of predicted sources if 
their SFR is statistically too high. Indeed, star formation 
processes critically depend on fundamental mechanisms 
and physical parameters that may still not be fully con- 
strained in the simulations, like for instance the critical 
density above which stars can form depending on the 
galaxy internal dynamics and the various energy feed- 
back (AGNs, supernovae, ...). 

Finally, it is worth noting that the observed optical and 
near-IR colors may also depend on the Initial Mass Func- 
tion (IMF) assumed in the simulations. In the model of 
Lacey et al. (2008) a normal neighborhood IMF leads 
to a too little evolution of the mid-IR luminosity func- 
tion compared to the data, and it also underestimates 
the strong evolution observed in the faint submillimeter 
number counts (Baugh et al. 2005). A key ingredient 
characterizing their scenario is the use of a top-heavy 
IMF for the bursts of star formation triggered by galaxy 
mergers, which seems to better reproduce the evolution 
of the galaxy luminosity function observed in the mid- 
IR. Under this assumption though, the larger fraction of 
blue and massive stars produced in bursts should lead 
to bluer stellar populations and this could be one of the 
reasons explaining the excess of blue galaxies in their 
predictions. A much more detailed comparison between 
the intrinsic properties (mass, SFR, SEDs) of individual 
simulated sources and those characterizing the observed 
24 jj,m galaxy population will be necessary to identify the 
critical parameters that need to be better constrained for 
improving model predictions in the infrared. 

Given the large number of ingredients to control, the 
determination of proper predictions of galaxy evolution 
at IR wavelengths with current semi-analytical models 
remains obviously a challenging task. Reproducing the 
evolution of the IR luminosity function requires not only 
an accurate description of the evolution of the bolomet- 
ric luminosity function for star-forming galaxies but also 
an accurate treatment of dust extinction as a function of 
galaxy properties and a reliable representation of galaxy 
SEDs in the infrared. Furthermore, dust-obscured star- 
forming galaxies are not the only contributers to the 
number counts at IR wavelengths. Quasars and dusty 
AGNs should also be taken into account, especially at 
the bright end of the luminosity function where their con- 
tribution is far from being negligible (e.g., Brand et al. 
2006). Given the coeval growth of bulges and super- 
massive black holes (e.g., Page et al. 2001) "composite" 
sources experiencing both nuclear and star-forming ac- 
tivity may actually represent a large fraction of the mid- 
IR selected galaxy population at high redshift. This large 
variety of physical properties characterizing IR-luminous 
galaxies added to the complexity of their spectral sig- 
natures in the mid-IR can thus make the predictions 
from semi- analytical simulations particularly difficult at 
IR wavelengths. 

This contribution of nuclear accretion to cosmic history 
has not been implemented in the model of Lacey et al. 
(2008). Therefore the feedback action of AGNs, which 
has been recently discussed in the framework of other 
semi-analytical models as an efficient process to regulate 
the star-forming activity in galaxies, is not taken into ac- 
count in their simulations. Instead, the decrease of gas 
cooling for the formation of massive spheroids is con- 
trolled in their model by the action of superwinds driven 
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by supernova explosions. These winds produce a feed- 
back effect qualitatively similar to that of AGNs, but the 
underlying physical mechanisms are somehow different 
and they result in different predictions for the evolution 
of the galaxy luminosity function with redshift (Lacey 
et al. 2008). Infrared simulations of distant sources ac- 
counting for the role of AGN feedback were actually pre- 
sented by Silva et al. (2005) using the model of Granato 
et al. (2004). However their scenario did not directly in- 
clude the merging of galaxies and dark matter haloes, 
and contrary to our results their 24 /xm source redshift 
distribution at S'24;^m > 0.2 mjy was predicted to peak 
at z^2. Given the current ACDM paradigm as well as 
the role that nuclear accretion played in the evolution 
of dusty high redshift galaxies, one possible approach to 
progress in our understanding of the deep IR Universe 
with semi-analytical models would be therefore to ex- 
plore the impact of scenarios including both the hierar- 
chical assembly of dark matter haloes and the contribu- 
tion of AGNs. 



8. SUMMARY 

We performed our very first analysis of the deep 24 y^m 
observations carried out with the Spitzer Space Telescope 
over the 2 deg^ of the COSMOS field. We detected al- 
most ~ 30 000 sources down to a 24 /xm flux of SOfiJy. 
Using the multi-wavelength ancillary data available in 
COSMOS we identified the optical/near-IR counterparts 
and the photometric redshifts for ^ 95% of these 24 /im 
detections. The main results that we derived can be sum- 
marized as follows: 

1 . The redshift distribution of the 24 /im selected 
galaxies shows a prominent peak at 2 ~ 1 but more than 
~ 50% of the sample is located at 1 < 2; ^ 3. Below 
z^l almost all 24 sources arc brighter than 24 mag 
(AB) in the i+-band, while a non-negligible of the MIPS 
sources at higher redshifts are associated with very faint 
optical counterparts (i.e., 24mag < z+< 26.5mag). 

2. Below a fixed optical magnitude limit the frac- 
tion of i+-band sources with 24 /zm detection strongly 
increases up to z^2. This reflects the strong evolu- 
tion that star-forming galaxies have undergone with look- 
back time. However this rising trend shows a clear 
break at z~1.3, which seems to be more pronounced 
when higher 24 /zm flux selections are considered (e.g., 
'S'24/im^0.3mJy). This feature is also visible in the red- 
shift distributions of 24/Ltm sources. It probably origi- 
nates from fc-correction effects implied by the presence 
of the strong aromatic features and/or the Silicate ab- 
sorption commonly observed in the spectra of luminous 
dusty galaxies. 

3. Using extrapolations of the number counts at faint 
fluxes we resolved the build-up of the mid-infrared back- 
ground across cosmic ages. We find that ~50% and 
~ 80% of the 24 fxm background intensity originate from 
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galaxies at z <1 and z <2 respectively, with a contribu- 
tion of faint som-ccs reaching ~20±10% below our de- 
tection limit. As already suggested by other previously- 
published studies we thus infer that the Cosmic Infrared 
Background is mostly produced by moderately-luminous 
dusty galaxies at 0.5<z<1.5 rather than very luminous 
and obscured sources at z>2. 

4. The comparison between our results and semi- 
analytical predictions reveal substantial discrepancies in 
the simulated redshift distributions as well as in the 
R — K colors of the 24 fim source counterparts. In partic- 
ular we found a strong excess of blue sources predicted 
at z ^ 2 and the K band fluxes of simulated galaxies ap- 
pear to be systematically too low at z>0.5. This could 
reflect an underestimate of the density of high redshift 
massive sources with strong on-going star formation pos- 
sibly related to fundamental physical parameters and/or 
processes that are still not fully constrained in the sim- 
ulations of galaxy formation. It could also be due to 
an inaccurate treatement of dust extinction in the mod- 
els, which would imply discrepancies in the predictions 
of the galaxy SEDs at rest-frame UV, optical and mid-IR 
wavelengths. 

5. Comparisons with some of the most recent back- 
ward evolution scenarios reproduce reasonably well the 
flux/redshift distribution of 24 /im sources up to z^3. 
Nonetheless none of them is able to accurately match 
our results at all redshifts over the full range of 24 /xm 
fluxes probed by our survey. 
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